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^ : ABSTRACT 

■ The magnetorotational instability (MRI) is the leading candidate for driving 

cri ' turbulence, angular momentum transport, and accretion in astrophysical disks. 

I consider the linear theory of the MRI in a thin, equatorial disk in the Kerr 

metric. I begin by analyzing a mechanical model for the MRI that consists of 

^ ' two point masses on nearly circular orbits connected by a spring. I then develop 

f^ ■ a local Cartesian coordinate system for thin, equatorial Kerr disks. In this local 

\Q • model general relativistic effects manifest themselves solely through changes in 

O 

O ■ MRI can be analyzed in the context of the local model using nonrelativistic 



(N 



^ ■ magnetohydrodynamics, and the growth rates agree with those found in the 

^^ mechanical model. The maximum growth rate measured by a circular orbit 



observer differs from a naive estimate using Newtonian gravity by a factor that 



O 

^ ■ varies between 1 and 4/3 for all radii and for all a/M. 



•^ ■ Subject headings: accretion, accretion disks, black hole physics, Magnetohydro- 

j_j ■ dynamics: MHD 



1. Introduction 

Balbus & Hawley (1991; hereafter BH) discovered a local instability of weakly magne- 
tized accretion disks, studied earlier in the context of magnetized Couette flow by Velikhov 
(1959) and others. This instability is widely known as the magnetorotational instability 
(MRI). The instability conditions (see Balbus & Hawley 1998) are easily satisfied, so almost 
every ionized astrophysical disk is likely subject to the instability. Numerical experiments 
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have demonstrated that the MRI can initiate and sustain magnet ohydro dynamic (MHD) 
turbulence in disks (e.g. Hawley, Gammie, & Balbus 1995). They have also shown that 
MHD turbulence in disks transports angular momentum outward. MRI initiated turbu- 
lence is therefore a plausible (arguably the most plausible) candidate for driving accretion 
in astrophysical disks. 

Disks around black holes are also likely subject to the MRI. The development of the 
MRI in these disks is of fundamental interest. But it is also of practical interest because 
numerical models of relativistic disks are now being developed (De Villiers & Hawley 2003; 
Gammie, McKinney, & Toth 2003; De Villiers, Hawley, & Krolik 2003; Gammie, Shapiro, & 
McKinney 2003). In particular, one would like to know if the maximum growth rate is very 
different from the (3/4)(GM/r^)^/^ expected from analyses of Newtonian disks. 

In this short paper I will study the MRI in thin, equatorial disks in the Kerr metric. I 
will use units such that GM = c = 1, where M is the black hole mass, and Boyer-Lindquist 
coordinates t, r, 9, (j), where 

rfs^ = - 1 - — dt^ + — dr^ + J: d9^ + — - — dc/)^ dcj) dt, (1) 

and Ti = r'^ + a'^ cos^ 9, A = r'^ — 2r + a^ and A = {r'^ + a^Y — a^A sin^ 9. It is useful to 
define the following relativistic correction factors that asymptote to 1 at large distance from 
the black hole: 

^ = 1 + aVr^ + 2aVr^ (2) 

i3 = l+ar-3/^ (3) 

C = 1 - 3/r + 2ar-3/^ (4) 

V = l-2/r + a^/r'^, (5) 

J^ =1- 2ar~^/^ + aVr^ (6) 

and 

^ = 1 - 2/r + ar-^/'^. (7) 

These definitions are identical to those used in Novikov & Thorne (1973). 

In what follows I will first develop a relativistic mechanical model for the MRI (§2), then 
a local Cartesian coordinate system that permits the treatment of the MRI near a rotating 
black hole using nonrelativistic MHD (§3). I compare the maximum growth of the instability 
to the shear rate in §4, and give a brief summary and prospects for future analysis in §5. 
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2. A Mechanical Analogy 

BH considered the linear theory of a magnetized disk in the WKB approximation. The 
full theory is rather complicated, but Balbus & Hawley (1992) developed a useful mechanical 
model for the instability that consists of two point masses on nearly circular orbit connected 
by a spring. The masses are analogous to fluid elements connected by a magnetic field, 
which behaves like a spring. This mechanical model captures key features of the full MHD 
instability, including the maximum growth rate and the fact that the instability grows by 
exchange of angular momentum. Indeed, when both the magnetic field and the wavevector 
of the perturbation are parallel to the axis of rotation the mechanical modes have frequencies 
which are identical to the WKB modes of the full MHD model if one replaces the natural 
frequency 7 of the spring-mass system by the Alfven frequency |k-VA|. It is therefore natural 
to begin a study of the MRI in the Kerr metric by considering a relativistic version of the 
Balbus & Hawley (1992) mechanical model. 

The orbits of free particles are described by the geodesic equation 

Uj iXj -t^II ^'^ HitXj / r\\ 

which has the following solution for circular equatorial orbits: 

«M = ^ = {BC-'/\ 0, 0, r-3/2c-^/2} (9) 

dr 

= M*{l,0,0,fi} (10) 

where Q = l/{r^/^ + a) = r-^/^B-\ 

Consider perturbations about the circular orbit, a;^ -^ x^ + ^^1 where ^^ is small. Then 



dr^ 



dj:iyu'C - 2r:;,tz^e' + o{m- (11) 



Substituting ^^ ~ exp(— icur) and evaluating the connection coefficients in Boyer-Lindquist 
coordinates, one finds 

oj\u^ - U^){0J^ - K^) = Q. (12) 

Here 



2 1 1 - Aa/r^/'^ + 3aVr2 



^ = r ^13) 



is the "vertical" frequency and 



. 1 1 - 6/r + 8ar-3/2 _ 3^2/^2 

" = 7i c <"> 
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is the epicyclic frequency. On the innermost stable circular orbit (ISCO) /t^ = 0. The other 
four modes are zero-frequency relabelings of the t and (f) coordinate. 

The frequencies u and k are defined with respect to the proper time r; they can be 
converted to frequencies with respect to the Boyer-Lindquist time coordinate t (and with 
respect to observers at infinity) by dividing by -u* = dt/dr = BC^^^"^ for a circular orbit. 
Using subscript t to denote a frequency with respect to t, 

nl = n^{l - 6/r + 8ar"^/^ - 3a^/r^) (15) 



-^ 



n\l-4ar-^/^ + 3ayr^). (16) 



Now consider two particles tied together by a spring. The particles lie at x'^ ± C,'^, and 
in the absence of rotation the natural frequency of the masses and spring is 7. The particles 
orbits now interact with the potential 

li'h,,ec (17) 

where /i^,^ = g^u + u^Uy projects the displacements into the space normal to the unperturbed 
(circular) orbit's four- velocity. 

The equations of motion can be derived in a variety of ways, perhaps most simply by 
varying the action with respect to ^^. The result is 



cir2 



d^T^^^u'^u'C - ^K.u'^e - YKC- (18) 



Assuming that S,^ ~ exp(— icjr), one finds that the nontrivial frequencies that correspond to 
motions in the equatorial plane obey 

where 

Equation (19) is quadratic in a;^ and so easily solved. By varying 7 one finds that the 
maximum growth rate is 

2 g^ _ 9 1 fvV 

"^max - 4(^2 + ^2) " J^^ I^^J ( ) 



which occurs for 

7^ = ^ V^- (22) 



s^ s2 + 2k^ 

4 S^ + K^ ' 
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The maximum growth rate is defined with respect to proper time r; it is the growth rate 
measured by an observer on a circular orbit at radius r. Converting to frequencies with 
respect to Boyer-Lindquist time t (observers at infinity): 

--L.. = l^' C^) (23) 



16 V C , 

This may be compared to the Newtonian maximum growth rate — c^maa; n — (9/16)r~^. I will 
refer to the ratio of the relativistic maximum growth rate to Newtonian maximum growth 
rate / = \uJmax/^max,N\ ^s the normalized growth rate. 

Figure 1 shows the run of the normalized growth rate for a = 0.9, 0, and —0.9 (i.e. a 
counterrotating disk). The most dramatic effect is the drop in normalized growth rate with 
respect to time t in the prograde disk. This is due to gravitational redshift (the conversion 
from time r to time t) so it is a global feature of the spacetime rather than a characteristic 
of the MRI: any disk instability would experience the same redshift. The effect is strongest 
on the ISCO for prograde disks. The normalized growth rate with respect to t drops to zero 
on the ISCO as a ^ 1, since -u* ^ oo there (see eq. 9). 

The very limited variation of normalized growth rate with respect to proper time r as 
r and a vary is peculiar to the MRI. This can be understood as follows. Since — u;^^^ = 
s^/(4(s2 + ^2)) and k^ = o at the ISCO, -cu^x = ^74 at the ISCO. One may also show 
that s^ + K^ = 4r~^, hence — u;^aa; = ''""'^ at the ISCO. Then tu^ax/'^max tv = 16/9 on the 
ISCO. Since the relativistic correction to the MRI growth rate is largest at the ISCO, the 
normalized growth rate varies between 4/3 at r = rjsco and 1 as r — > oo. 



3. Local Model 

In studies of thin disks in a Newtonian potential it has proven useful to introduce a local 
Cartesian coordinate frame, or local model. This frame is centered on a circular orbit, the 
X axis points along the radius vector, the y axis points forward in azimuth, and the z axis 
points up, normal to the disk. Beginning with a global spherical coordinate system r, 9, 0, 
the transformation to nonrelativistic local coordinates is 

X = r — tq (24) 

l/ = ro(0-0o-fit) (25) 

where 0o is a constant and here only Vt = r^ , and 

z = ro(|-^). (26) 
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Expanding to lowest order in |x|/ro ~ |x|/rof2, the free particle equations of motion are 

x = 2VLy + Sfi^x, (27) 

y = -2nx, (28) 

i = -n^z. (29) 

The linear theory of magnetized disks in the WKB approximation can be carried over in its 
entirety to this local model because higher order terms in e are also higher order terms in 
the WKB approximation. 

In a relativistic disk in the equatorial plane of the Kerr metric we may proceed in 
exactly the same way. We want to erect local Cartesian coordinates near a circular orbit. 
In Boyer-Lindquist coordinates the circular orbit is located at r = tq and has d(j)/dt = VLq = 
l/(rQ + a). The transformation to local coordinates r, x, y, z is 

r = tgC-^l^ - (0 - 0o)ro'/'^C-^/2 (30) 

(31) 

(32) 

(33) 

This transformation is determined by the conditions that the local coordinates should be 
locally Minkowski, and that the point x,y,z = should follow a circular orbit. This trans- 
formation is (not by accident) identical to the transformation to a corotating orthonormal 
basis of one- forms given by Novikov & Thorne (1973). 

Now consider the limit e^ = {x^ + y"^ + z^)/i^o ^ 1- Expanding gu to O(e^), gu to (9(e) 
and gij to 0{1) (higher order terms are not required for a consistent expansion of ds"^; one 
can think of dr as being 0{1) and rfx* as being 0{e)), the metric becomes 

ds"^ = (-1 + s^x^ - iy^z'^)dT^ + AVtxdrdy + dx^ + dy"^ + dz^ (34) 

where Vt = r^ , 

V 



X = 


(r- 


- ro)V-'/', 




ro{(p- 


-00- 


- Qot)BV^/^r 


-1/2 


z 


= ^0 


il-er 





and 



^ = ^[J:\' (35) 

.= = -Lfi^i<!^±Mi|. (36) 



These are the same as the s and v defined in the last section. 
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The metric component Qrr can be interpreted as ^ —1 + 2-?/' where ip is the effective 
gravitational potential, and the Qry term indicates that the local model is a frame rotating 
with frequency Cl. Using this interpretation, or working directly with the geodesic equation, 
one arrives at the equations of motion for a free particle 

x = 2Qy + s^x, (37) 

y = -2nx, (38) 

z = -u^z, (39) 

where we have assumed that (x^ +y'^ + P)/{QrQ) ~ e^. One can look for motions with time 
dependence exp{—iujT); the free particle frequencies are identical to those found from the 
perturbed geodesic equation in §2. 

If the disk is thin then the sound speed is small compared to the speed of light. Since one 
expects the fluid speeds within the disk to be subsonic, the fluid motions are nonrelativistic. 
It follows that in the local Cartesian frame the fluid obeys the equations of nonrelativistic 
fluid dynamics. The MHD equations of motion in the local model are then 

Dv ~ 9 9 Vp V52 (B-V)B , , 
= -2fl X V + s^xe^ - u^ze^ h ^ —, (40) 

where fl = Qsz. These equations, coupled to the continuity equation, the (nonrelativistic) 
induction equation, an energy equation, and the condition V ■ B = 0, give a complete 
description of the local evolution of the plasma in the MHD approximation. 

The nonrelativistic linear analysis of the Balbus-Hawley instability can now be carried 
over into the local frame by replacing the Newtonian rotation frequency Q with Cl, the 
coefficient of the tidal expansion of the effective potential (3/2)f2^ by s^, and the Newtonian 
vertical frequency Q by z/. General relativistic effects appear solely through the Coriolis 
parameter Cl and the tidal expansion parameter s. 

To carry forward a linear analysis we need to specify an equilibrium. The simplest choice, 
which captures most of the features of the nonrelativistic instability, is a uniform vertical 
magnetic field near the midplane of the disk, so the Alfven velocity can be written v^ = 
(i?/(47rp)^/^)e^. I will restrict consideration to perturbations of the form exp(— icur + ikzz) 
(more general perturbations do not introduce qualitatively new features in ideal MHD). Then 
the dispersion relation is 

^4 _ ^2(^2 ^ 2(k ■ v^)2) + (k ■ v^)2((k ■ v^)^ -s')=0 (41) 



Since this is identical to the relation obeyed by the mechanical model, with 7^ ^ (k ■ ^aY-, 
the maximum growth rate of the MHD instability is the same as that of the mechanical 
model. 



4. Relationship to Shear Rate 

Balbus & Hawley (1992) conjectured that the maximum growth rate should always be 
half the shear rate. To check this, we need a relativistic generalization of the shear rate. The 
natural generalization involves the scalar 



(T^ = Or^'''(J^^u 


(42) 


where 


(43) 


is the rate-of-strain tensor and 


(44) 


is the divergence of the velocity field. For a Cartesian shear flow Vy ■ 


= Ax, where A is 



constant, one finds o"^ = A^/2, so 2(7^ is the square of the shear rate. 

For circular equatorial geodesies in the Kerr metric a^ can be evaluated in any frame, 
including the Boyer-Lindquist coordinate frame (this is a lengthy calculation!). Novikov 
& Thorne (1973) evaluated a^i, in an orthonormal tetrad moving with the flow, which is 
equivalent to our local coordinate frame. They found 



(y{r){4>) = ^W(r) = ^^ ( -^ J , (45) 

and all other components are zero. Then 

Evidently 

-^L. = ^(2^') (47) 

This confirms the relationship between the maximum growth rate and the shear rate sug- 
gested by Balbus & Hawley (1992). In the local coordinate system, of course, the relationship 
is trivially verified. 
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5. Conclusion 

I have shown that the growth rate of the MRI in thin, equatorial disks around Kerr 
black holes does not differ sharply from its Newtonian counterpart, at least as measured by 
an observer on a circular orbit. As measured by an observer at infinity the growth rate is 
redshifted strongly for instabilities growing near the marginally stable orbit in corotating 
disks. 

A similar analysis could be carried out for inclined disks. The work of Bardeen & Pet- 
terson (1975) and successors (see Demianski & Ivanov 1997 and references therein) suggests, 
however, that close to the hole disks settle into the equatorial plane of the black hole in 
the presence of angular momentum diffusion (such as that which might result from MHD 
turbulence initiated by the MRI). At larger distance the deviations from Newtonian gravity 
will be small and will likely not affect the development of the most unstable modes, which 
grow on the dynamical timescale. 

Earlier work by Gammie & Popham (1998) (see footnote 7) stated the results presented 
here without proof. Araya-Gochez (2002) gives an impressive fully relativistic linear analysis 
that is equivalent to results presented here (compare Figure 1 with his Figure 2). The fully 
relativistic analysis is perhaps the ultimate proof that the MRI carries over nearly unchanged 
to relativistic disks; the present paper states these results in a greatly simplified form. 

This work was supported by NSF ITR grant PHY 02-05155 and by NSF PECASE 
grant AST 00-93091. I wish to thank John Hawley and an anonymous referee for helpful 
comments. 
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Fig. 1. — The figure shows the maximum growth rate of the Balbus-Hawley instabihty for 
a thin disk in the equatorial plane of the Kerr metric normalized to the naively calculated 
nonrelativistic growth rate. The left panel shows /,-, the normalized growth rate with respect 
to the proper time r of an observer on a circular orbit. The right panel shows ft, the 
normalized growth rate with respect to the Boyer-Lindquist time t, which is the growth rate 
measured by an observer at large radius. 



